Finite Temperature Phase Diagram of Quasi-Two-Dimensional Imbalanced Fermi 

Gases Beyond Mean-Field 
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We investigate the superfluid transition temperature of quasi-two-dimensional imbalanced Fermi 
gases beyond the mean-field approximation, through the second-order (or induced) interaction ef- 
fects. For a balanced Fermi system the transition temperature is suppressed by a factor « 2.72. For 
imbalanced Fermi systems, the polarization and transition temperature of the tricritical point are 
significantly reduced as the two-body binding energy |es| increases. 
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I. INTRODUCTION 

The domain of advanced and accurate experimental 
techniques of laser cooling and magnetic trapping, as 
well as imaging of neutral ultracold two-spin-component 
atomic Fermi gases, has permitted several crucial in- 
vestigations in these systems. One of the very impor- 
tant studies achieved is the crossover from the Bardeen- 
Cooper-Schrieffer (BCS) phase of weakly bound Cooper 
pairs to the Bose-Einstein condensate (BEC) phase of 
strongly bound diatomic molecules in three-dimensional 
(3D) trapped Fermi gases 

The lowering of dimensionality in many-body quantum 
systems opens a window for the appearance of interest- 
ing phenomena, such as the (yet unexplained) effect of 
high-temperature superconductivity in two-dimensional 
(2D) cuprates. Thus, the trap geometry, which is also 
currently under full control, is fundamental in the use 
of ultracold atomic gases in the simulation of condensed 
matter systems, such as layered 2D strongly correlated 
superconductors. 

As examples of very recent experimental achievements 
in cold balanced Fermi systems in low dimensions, more 
specifically in 2D, we highlight the report of the measure- 
ment of the density profile and temperature of a 2D gas 
of atoms [i[ , pairing in a harmonically trapped 2D atomic 
Fermi gas in the regime of strong coupling [5[ , and the de- 
tection of a many-body pairing gap above the superfluid 
transition temperature (the pseudogap phenomenon) Q . 
To the best of our knowledge, up to now the only exper- 
imental investigations with imbalanced 2D Fermi gases 
concerns the Fermi polaron problem, in which a single 
spin-j. atom interacts strongly with a Fermi sea of spin-'f 
atoms 0- 

According to the Mcrmin-Wagncr-Hohcnberg- 
Coleman (MWHC) theorem Q, the long-range order is 
destroyed by fluctuations in uniform, 2D systems. This 
prohibits the formation of a superconducting phase with 
a homogeneous order parameter, associated with the 
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breaking of a continuous symmetry. Nevertheless, 2D 
systems may undergo a phase transition to a state with 
quasi-long range order via the Berezinskii-Kostcrlitz- 
Thouless (BKT) transition Q. Since this transition docs 
not require symmetry breaking, it is not forbidden by 
the MWHC theorem. 

By taking into account the phase fluctuation effect, the 
BKT transition has been investigated in 2D balanced [13, 
ITiT ] as well as imbalanced systems [l2j ■ In all these works 
a strong dependence of the BKT transition temperature 
Tbkt on the two-body energy |e#| has been found. Below 
?bkt pairs of vortices and antivortices surge, and will 
eventually condense as the temperature is lowered fiol - 

The direct consequence of the MWHC theorem is that 
one should find T c = for the phase transition tempera- 
ture between the normal and superfluid phases in a pure 
2D system. However, quasi-2D systems (as, for exam- 
ple, a stack of planes where tunneling between different 
planes is completely suppressed by a large trapping po- 
tential, constituting effective 3D systems) may exhibit a 
superfluid phase transition at a finite T c with no conflict 
with the MWHC theorem [II]. 

In this paper we study the finite-temperature (T) 
ground state of a quasi-2D atomic Fermi gas with chem- 
ical potential imbalance, beyond the mean-field approx- 
imation. The standard mean-field calculation does not 
take into account the effects of the medium on the two- 
body interaction. This correction, considered first by 
Gorkov and Mclik-Barkhudarov (GMB) [14(, has been 
referred to as induced interactions [HI , and was found to 
suppress the mean-field critical temperature of a 3D bal- 
anced Fermi gas by a factor sa 2.22 [T3 - [l7t . In addition, 
it has been shown that the GMB correction substant ially 
reduces the order parameter in 2D and 3D lattices [19| . 
The influence of induced interactions in a 3D imbalanced 
Fermi gas has been taken into account in Rcf. [2(|, and 
it was found that the polarization P and the transition 
temperature T of the tricritical point are both reduced 
from the mean-field results by a factor of about 2.22, 
meaning that the transition temperature suppression is 
comparable with the one found for the spin-balanced con- 
figuration. 
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The possibility of the Fuldc-Fcrrel-Larkin-Ovchinnikov 
(FFLO) state with modulated order parameter [2l[ is 
ignored in this work. As in [T34l9l |. we consider only 
pairing between atoms with equal and opposite momenta. 

The paper is organized as follows: In Sec. [TT| we ob- 
tain the finite temperature thermodynamic potential of 
the model in the mean-field approximation. In Sec. IIIII 
we review some basic zero temperature mean-field results 
which will be used latter. In Sec. II VI we construct the fi- 
nite temperature phase diagram beyond the mean-field 
approximation through the GMB correction. We con- 
clude in Sec. |Vj 



II. THE MODEL HAMILTONIAN AND THE 
MEAN-FIELD THERMODYNAMIC POTENTIAL 

We start by considering a 2D nonrelativistic dilute sys- 
tem of fermionic atoms of mass M, with two hypcrfine 
states labeled as a =t,4~ This spin-t and -\. mixture 
could, in principle, be obtained in 2D experiments with 
the two lowest hypcrfine states of 40 K la, @], or with 6 Li 
atoms, as in the 3D experiments [22H24| . Their single- 
particle dispersion relations are given by = 2M . 
Throughout the paper we set h = 1. The quasi-2D sys- 
tem can be modeled by the following pairing Hamilto- 
nian: 



H = H„ 



where 



H = J2 4ri(k)Mk), 



(1) 



(2) 



=t4 



is the kinetic (free) part of H and Hi nt is given below. 
ipl{k) and Vv(fc) in Eq. © are the creation and annihi- 
lation operators, respectively, for the f and -J, particles. 
To assure population imbalance in the system, we have 
introduced different chemical potential for the species a 
as jU CT = /i + ah, where ah = ±h. Then, the chemical 
potentials fi a fixes the number densities n a of the differ- 
ent fermions. The new dispersions for the free species a, 
relative to their Fermi energies, are e£ = — fj, a . The 
interaction Hamiltonian is given by 



Hi, 



gY,^)i>\{-k)^i{-k')^{k'), (3) 



k.k' 



where the bare coupling constant g is negative, to express 
the attractive (s-wave) interaction between the spin-'f 
and -J, fermionic atoms. 

After the mean-field (MF) approximation and the di- 
agonalization of the expression for Hmf, we arrive at 
the following expression for the thermodynamic poten- 
tial, from which all thcrmodynamical quantities of inter- 
est can be obtained: 



n = 



k 



E k - Tln(e 



-/3£? 



Tln(e 



(4) 

1) 



where, for simplicity of notation we have labeled 4-= a, 
f = b. Here (3 = 1/T, where we have set the Boltzmann 
constant equal to 1. We have also defined £^' b = Ek±e k 



as the quasiparticle excitations, with Ek = y + A 2 , 

e k = f k 2 k ' an< ^ ^ ne constant pairing gap is given by 
A = -5/(0(^(-^(fc))=A*. 

To regulate the ultraviolet divergence associated with 
the zero temperature term in Eq. (|4]), we introduce the 
2D bound-state equation (2j| 



1 
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d 2 k 



(2tt) 2 2^ + \e B \' 



(5) 



where es is the 2D two-body binding energy. In addition 
to the regularization of the ultraviolet divergence present 
in Eq. ((4]), and consequently in the gap equation (sec 
below in Scc. UVBj) . this equation relates the strength g of 
the contact interaction with e b , which is more physically 
relevant, as will be clear now. In order to make contact 
with current experiments, it is convenient to relate cb 
to the three dimensional scattering length a s . In the 
scattering of atoms confined in the axial direction by a 
harmonic potential with characteristic frequency wt they 
are related by [H, [2?| 



Cuj, 



-exp 



Z7T 

a s 



(6) 



where a s is the 3D s-wave scattering length, u>l = 
y/8n 2 V /(M\ 2 ), l L = 1/y/MJE, and C w 0.915. V Q 
is the amplitude of the periodic potential Vo sin 2 (27rz/A) 
generated by two counter-propagating laser beans with 
length A parallel to the z axis [26[ . 



III. BASIC MEAN-FIELD ZERO 
TEMPERATURE RESULTS OF 2D 
IMBALANCED FERMI SYSTEMS 



In this section we borrow from Refs. [25J, |28|, |29| some 
basic zero-temperature MF results, which will be needed 
in the next section. 



A. Balanced systems 

The gap and number equations are obtained by 
dfl/dA = and n a = —dtt/dfj, a , respectively. For the 
balanced system where h = 0, we find 
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s/WT&-li = \e B \, (7) 

and 

v / M 2 + A 2 + A i = 2e F , (8) 

where the two-dimensional Fermi energy is defined as 
ep = ^jj-, with riT = n a + rifa. In the balanced configura- 
tion n a = ni, = n = ^f- . Solving these two equations self- 
consistently we arrive at the well-known results [25|, H^] 



A = y/2e F \e B \, (9) 

and 



where we are considering fi a and h positive. 

From the graphical analysis of f2 as a function of A for 
various asymmetries [H, [30, [H| , one sees that increasing 
the imbalance h and keeping \i fixed, the minimum is still 
located at Ao up to a maximum or critical imbalance 
h c , after which there is a quantum phase transition to 
the normal state with A = 0. h c is found through the 
equality Hq = 0,^ nb , which yields 

fcc = Mo|e B | + (y) 2 . (15) 

From this one easily finds, by plugging into Eq. [15] fio 
from Eq. in the BCS limit (|e B | « e F ): h c = 
which is the same 3D result known as the Chandrasekhar- 
Clogston limit of superfluidity [lU, [33| . 



Mo = £f r - . (10) 

The value of the free energy at the minimum is 



n(h = o,A = A ) = n = - K \^ + ^j , (ii) 

where k = y~ , whereas the energy of the balanced normal 
state is given by 

Cl(h = A = 0) = QZi = -«Mo- (12) 

A direct comparison between Eqs. (fTTj) and (|T2|) shows 
that the superfluid state is energetically preferable to the 
normal state for any eg ^ 0. Since a two-body bound 
state exists even for an arbitrarily small attraction in 
2D [2|| , the pairing instability will always happens in 2D 
balanced two-component Fermi systems at T = 0. 

B. Imbalanced systems 

We now turn our attention to the cases where h 7^ 0. 
The free energy of the unbalanced normal state. Q(/i, A = 

°) = n iL,: is s ivcn b y 

f k F d 2 k /" fc ^ d 2 k 

Qimb= J j^y tl+ L w? el (13) 

which gives, after the integration in k, the free energy 
of a (normal) two species gas of fermionic atoms in two 
dimensions 



IV. THE PHASE DIAGRAM BEYOND THE 
MEAN FIELD 

Now we construct the phase diagram of the model 
at finite temperature considering corrections beyond the 
mean field, taking into account the GMB correction. 

A. Induced interaction in a spin polarized Fermi 

gas 

The induced interaction was obtained originally by 
GMB in the BCS limit by the second-order perturbation 
[l5| . For a scattering process with P1+P2 — > P3 +P4, the 
induced interaction for the diagram in Fig. 1 is expressed 
as 

Uind(pi,P4) = -g 2 X P h(pi - Pi), (16) 

where pi = (k j , ui i ) is a vector in the space of wave- vector 
k and fermion Matsubara frequency w; = (21 + 1)-k/(0). 
Including the induced interaction, the effective pairing 
interaction between atoms with different spins is given 
by 




FIG. 1: The lowest-order diagram representing the in- 
duced interaction f7i n d(pi,P4). Arrowed and dashed lines de- 
scribe fermionic propagators and the coupling g between the 
fermionic atoms. 



= -^1 + ^1], (14) U ef f(p u P4.) = U eff = g + U iaA (p uPA ). (17) 
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The polarization function Xph(j>') is given by 



of an imbalanccd Fermi s function of its spin po- 

larization P: 



X P h{p') 



h 2 /3A 



^2 50b{p)90a(p + P) 



d 2 k /] 



fS 

J k+q 



(18) 



where p' = (q,Qi), flj = 21tt//3 is the Matsubara fre- 
quency of a boson, A is the area of the system, and f(k) is 



/■it:;' 



1 



the Fermi distribution function /(£?' ) = l/(e 
The Matsubara Green's function of a non-interacting 
Fermi gas is given by Soo-(p) = — eko-)- The static 

polarization function is then 



%,/.(<?) = (0) 




(22) 



where 7 = 1 + cos <j> and 9p = 0(cos <fi — vl — -P 2 )- 

The s-wave part of the effective interaction is approx- 
imated by averaging the polarization functio n Yvh( q), 
which means an average of the angle <fi [3, EH Il7l |2C| : 



(x P h(q)) 



2tt 



d(f> X P h{q) = X P h(h). (23) 



„ f d 2 k 

X P h(\q\) = -2m J 



fa 
Ik 



(2-7I-) 2 q 2 - Amh + 2kq cos 6 
■b 



(19) 



fk 



q 2 + Amh — 2kq cos ( 



where q = \q\. At zero temperature the real part of 
Eq. (fT9|) is given by 



X P h(q) 



(20) 



-N(0) 



-N(0) 





4mhy (2k b F 



for q>k a F + k b f 



= -N(0),foTq<k F + k F , 
is the 2D density of states at the Fermi 



where N(0) 
level. 

In the scattering process the conservation of total mo- 
mentum implies that k\ + = k% + k^, with ki = —k% 
and k 3 = — k±. q is equal to the magnitude of k\ + £3, 



so q 



{k 1 + k 3 ).(k 1 +k 3 ) 



k\ 



k 2 



2k x .k z 



fcf + 2|fci||fc 3 | cos0, where <j> is the angle between 



ki and /C3. Since both particles are at the Fermi surface, 
\h \ = \k 3 \ = k F = y/2Mji, thus, q = fc Fv /2(l + cos 4>). 



The spin polarization is defined as P = "^ + " a 
from Eq. (| 14[) we find n Q = (M/27r)^ Q , yielding 



P 



/? 

/'' 



Then, 



(21) 



We are considering < n a < n^, which gives < P < 1. 
This is translated into having < h < fx. Equation (12"T]) 
will enable us to express the polarization function Xph(h) 



B. The GMB correction to the mean-field 
transition temperature and tricritical point 

Now we use the induced interaction effects calculated 
in Sec. IIV Al in order to obtain the corrected beyond- 
mean-field transition temperature of a balanced Fermi 
gas, and the tricritical point, Ptc = (ht c ,Tt c ), of an un- 
balanced Fermi gas. 

Extremizing the grand potential of Eq. (j4]) with respect 
to A and passing to integrals, we obtain a gap equation. 
The critical temperature T c is, by definition, the temper- 
ature at which A = 0. Then we find: 



- - Xph(h) 
d 2 k 1 



(2tt) 2 2e k 
d 2 k 1 



(24) 



(2tt) 2 2^ k 
d 2 k 1 
(2tt) 2 2~T k 



[l-jK)-/(4)] 
Xphih) 



[1 - M) f(4)] = 0, 



where e k = k 2 /2m — \i. We have used Eq. ([5]) to regu- 
late the ultraviolet divergence, associated with the zero 
temperature logarithmically divergent term in Eq. (|24p . 
The difference between the equation above and the usual 
mean-field thermal gap equation is that the particle-hole 
fluctuation has been taken into account through the ef- 
fective s-wave interaction U e ff. 

Thus, to find the transition temperature of a balanced 
Fermi gas below which the formation of Cooper pairs be- 
comes favorable, we have gone beyond the simple BCS 
approach and considered the second-order interaction ef- 
fects calculated in Eq. ([20]) . The fact that both scatter- 
ing particles are at the Fermi surface restricts the angle 
between k\ and k 3 to the single value <j> = 0, yielding 
q = 2kp. This gives Xph{h = 0) = —N(0), which results 
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T C ,GMB 



T C ,MF ^ T C ^ IF 
e ~ 2.72 : 



(25) 



agreeing completely with a previous work performed in a 
slightly different manner, but also considering the GMB 
correction [l^- The result above shows that the sup- 
pression of the transition temperature by second order 
interaction effects is higher in 2D than in 3D, where 

t^mb « iff EM3- 

In oder to obtain the GMB correction to the mean-field 
tricritical point, we define from Eq. 



f(T,h) 



d 2 k 



1 



(2tt) 2 2^ + \e B 
d 2 k 1 



- X P h(h) 



(26) 



(2tt) 2 4e fc 



tanh(^) 



tanhfe 



A graphical inspection shows that at the (tri) critical 
chemical potential imbalance h tc , both f(T, h) and its 
derivative with respect to T, g(T, h) = df(T, h)/dT, are 
zero [3~i| . The vanishing of these two functions corre- 
sponds to a = f3 = 0. a and j3 are the first and sec- 
ond coefficients of the expansion of the free energy in 
terms of the gap parameter, according to the Landau 
theory of phase transitions. By this criterion, one finds 
h tc and T tc of the tricritical point. Below this point, 
i.e., at low temperatures, the transition is of first or- 
der and the critical temperature has to be found by 
properly equating the energies of the normal and super- 
fluid phases. Thus, for a given h tc < h < h c , where 
h c is the Chandrasekhar-Clogston limit of superfluidity 
mentioned just below Eq. (|15|) . the temperature where 
n(A = 0,T) = Q(A (T),T) corresponds to the first or- 
der phase transition critical temperature. 

Our results for the induced interaction corrections to 
the mean-field tricritical polarization and temperature 
are shown in Figs. [2] and [3l respectively. As one can 
see in these figures, the second-order interaction effect 
in the BCS regime reduces the tricritical point by ap- 
proximately the same amount found analytically for the 
balanced scenario. It is easy to verify that in the BEC 
region, the effect of the induced interaction vanishes due 
to disappearance of the Fermi surface. From Eq. (|20l) we 
obtain that Xph(q) goes to zero in the limit k F = k b F = 0. 



V. CONCLUSIONS 

We have studied the effects of the induced interactions 
on the transition temperature of a quasi-2D imbalanced 
Fermi gas of atoms in the BCS region. We find that the 
transition temperature is reduced by a factor ps 2.72, in 
the case of a balanced Fermi system, and the temperature 
and polarization of the tricritical point of imbalanced 
Fermi systems, are also suppressed by approximately the 
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FIG. 2: The critical polarization P c = — of the tricritical 
point plotted as a function of the two-body binding energy. 
The solid line shows the MF result, and the dashed line shows 
the results corrected by second-order interaction effects. 
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FIG. 3: The tricritical temperature Tt c plotted as a function 
of |es|//^. The solid line shows the MF result, and the dashed 
line shows its corrections by second-order interaction effects. 



same factor as a function of |es|//Lt, in comparison to the 
MF results. 

It is worth mentioning that in the BEC region the uni- 
versal behavior found for the transition temperature in 
2D fermionic systems, considering the phase fluctuation 
effects [Tol - fl^l ]. is that the transition temperature obeys 
a behavior which may be expressed qualitatively in the 
form T c jtp rj const x tanh(|eB|/ep), reaching a limiting 
value for large \cb\/^f- The value of the constant is 0.075 
in Ref. [HI, and 0.125 in Refs. [rJ and [H]. Monte Carlo 
simulations give T c /ef ~ 0.1 for \cb\/^f = 10 [35] ■ This 
behavior can be explained with very simple arguments. 
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The phase fluctuation treatment correctly captures the 
physics that in the BEC side (large values of \es\) the 
system behaves as weakly interacting dimcrs of fcrmionic 
atoms, such that the transition temperature turns out to 
be very insensitive to the intensity of the interaction be- 
tween these composite bosons. 

We remark that to study the whole BCS-BEC 
crossover in a strictly 2D system, more detailed calcu- 
lations are required to take into account the effects of 
both second order interactions and phase fluctuations. 

In conclusion, taking into account the induced inter- 
actions, we have found the tricritical temperature and 
polarization of the superfluid phase transition of two- 



component quasi-2D imbalanced Fermi gases. Our re- 
sults are promising for achieving this transition in the 
regime of BCS pairing in current experiments. 
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